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Abstract. We use a simple tight-binding (TB) model to study electronic properties of free graphene flakes. 
Valence electrons of triangular graphene flakes show a shell and supershell structure which follows an 
analytical expression derived from the solution of the wave equation for triangular cavity. However, the 
solution has different selection rules for triangles with armchair and zigzag edges, and roughly 40000 atoms 
are needed to see clearly the first supershell oscillation. In the case of spherical flakes, the edge states of 
the zigzag regions dominate the shell structure which is thus sensitive to the flake diameter and center. A 
potential well that is made with external gates cannot have true bound states in graphene due to the zero 
energy band gap. However, it can cause strong resonances in the conduction band. 

PACS. 73. 21. La Quantum dots - 81.05.Uw Carbon, diamond, graphite - 61.48.De Structure of car- 
bon nanotubes, boron nanotubes, and closely related graphitelike systems - 81.05.Uw Carbon, diamond, 
graphite 



1 Introduction 

Electrons confined in a finite cluster of atoms with spheri- 
cal symmetry exhibit a shell structure |ll2j . In large enough 
clusters, the shells representing different classical periodic 
orbits can interfere forming a supershell structure [3 4 that 
has been observed in large alkali metal clusters |;5j. The su- 
pershell structure is especially visible in a two-dimensional 
triangular cavity which has only two classical periodic 
orbits^. The triangular cavity is interesting also due to 
the fact that the Schrodinger equation and the wave equa- 
tion are exactly solvable in that svstem [7l8l9j . and it has 
been shown that triangular shapes are preferred in two- 
dimensional nearly free electron svstems |10|ll|12| . 

Recently, experiments have shown that single layer 
graphene flakes can be prepared on inert surfaces where 
the graphene-surface interaction is weak fl 3ll4ll5ll6ll7llj^ 
Since the manipulation of graphene on different substrates 
is still a fast developing area, it is not out of question that 
graphene flakes with accurate shape and size can be even- 
tually processed on a substrate where the interaction is so 
weak that it does not affect the graphene electronic lev- 
els close to the Fermi point. Hence, we study ideal free 
graphene flakes neglecting the interaction with the sub- 
strate. The experiments have inspired a wealth of theoreti- 
cal studies of graphene [19 20 21 22 23 24 25 26 27 28 29 30 
but ackording to our knowledge the shell and supershell 
structure of large graphene flakes has not been addressed 
except in our recent work [3.3j . 

Electronic structure calculations based on the density 
functional theory (DFT) have shown that the energy levels 
close to the Fermi level, which consists of discrete points in 
graphene, are determined by the p electrons perpendicular 



to the graphene plane (for a review see [34]). A simple 
tight-binding (TB) model with only one electron per site 
and only the nearest-neighbour hopping describes well the 
electronic structure close to the Fermi points as suggested 
by Wallace already in 1947 [35] ■ The TB hamiltonian used 
is then the simple Hiickel model 

_ f — i, if i,j nearest neighbours 

~ \ 0, otherwise, ^ ' 

where the hopping parameter t (resonance integral) de- 
termines the width of the bands and the on-site energy is 
chosen to be = 0. We present our results in units t — 1 
which in real graphene corresponds to ~2.6 eV. In reality, 
the flake edges are either passivated (e.g. with hydrogen) 
or reconstructed in order to remove dangling bonds. The 
passivation is not expected to affect the perpendicular Te- 
states, and we can simply neglect the existence of such 
atoms. This has been also validated by our recent DFT 
calculations [33]. At the bottom of the valence band the 
TB model results in free-electron-like states with nearly 
constant density of states (DOS). This allows us to com- 
pare these "normal" free-electron states with those of the 
"massless electrons" at the bottom of the conduction band 
3Pgif ulated for the exactly same geometry. 

The paper is organized as follows: In Section 2 we 
discuss the shell and supershell structure in triangular 
graphene flakes, in Section 3 we show how the edge geom- 
etry dominates the shell structure in circular flakes, and in 
Section 4 we describe quantum dots that have been made 
with external potentials in an inflnite graphene sheet. Sec- 
tion 5 gives the conclusions. 
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2 Shell structure of triangular graphene flakes 

In a triangular cavity with hard walls and a constant po- 
tential inside, the two-dimensional Schrodinger equation 
has an exact solution [B] with energy levels 



/ 2 



(2) 



where eo depends on the particle mass and the size of the 
cavity, and m and n are integers with n > 2m > 1. Fig- 
ure [T] shows the density of states calculated with the TB 
model (TB-DOS) for a large graphene triangle of 44097 
atoms which has a zigzag edge. In such a large triangle, 
DOS is similar to that of an infinite graphene sheet ex- 
cept for the appearance of the edge states which appear 
as a sharp peak at zero energy. A detailed study of the 
energy levels at the bottom of the valence band reveals 
that the level structure is nearly exactly described with 
the analytical formula of Eq. TB-DOS and Eq. ^ 
produce the same curve shown in the lower left corner of 
the figure. Note that TB-DOS is plotted here as a func- 
tion of the wave number defined as Q = ^/e + 3t (The 
bottom of the band is — 3i). The regular oscillation as a 
function of Q corresponds to the shell structure and the 
peak amplitude variation (breathing) marks the supershell 
structure [6 . Triangles with an armchair edge show a sim- 
ilar supershell structure at the bottom of the band|33j. 




Fig. 1. Upper panel: TB-DOS (pz electrons) of a graphene tri- 
angle with 44097 atoms and zigzag edges. The discrete energy 
levels have been smoothened with Gaussians. The peak at zero 
energy corresponds to the edge states. The lower panels show 
TB-DOS as a function of the wave number at the bottom of 
the band (left) and above the Fermi level (right). The dashed 
line shows the analytical result of Eq. ([sjl . The wave numbers 
are: Q = y/e + 3 and q = e, where e is the energy in units of 
t = 1. The Gaussian widths have been adjusted in the lower 
panels in order to show the individual energy levels. 



At the bottom of the conduction band, close to e = 0, 
the dispersion relation of the electron energy is linear, 
e(fc) = cfik, where c is the electron velocity. This means 



that the electrons behave as massless particles. If we con- 
sider the conduction electrons as free particles, we cannot 
solve the energy eigenvalues from the Schrodinger equa- 
tion but should use the relativistic Dirac equation [34] . 
However, we choose here a simpler approach and make 
an ansatz that the energy eigenvalues are solutions of the 
Klein-Gordon wave equation with positive energy eigen- 
values. This immediately gives 



ei 



(3) 



which is the same as for normal electrons apart of the 
square root dependence of the quantum numbers and a dif- 
ferent prefactor. The numerical solutions of the TB prob- 
lem for triangles with a zigzag edge, indeed, show that 
the energy eigenvalues become more-and-more accurately 
described with those of Eq. ^ when the triangle size in- 
creases. Figure fl] (lower right panel) compares the result 
of the ansatz of Eq. ([3| with the full TB calculation for 
the large triangle. The agreement is nearly perfect up to 
e — 0.1 (corresponding to 0.25 eV), and the discrepancy 
at larger energies is due to the increasing nonlinearity and 
anisotropy of the energy bands. 

The results suggests that the supershell structure of 
the triangular cavity appears in zigzag-edged triangular 
graphene flakes, but the flake should have at least ca. 
40000 atoms (i.e., L > 40 nm) before the first supershell 
oscillation becomes clearly visible. We remark that in the 
case of armchair edge, Eq. (|3l is still valid, but also in- 
dices with m = n are allowed [55]. More detailed results 
for smaller triangles are described in Ref. [5H], where we 
also show that the shell structure is quite robust against 
edge roughness in the close vicinity of the Fermi level. 



3 Shell structure in circular graphene flakes 

In the case of a two-dimensional cavity with circular sym- 
metry, the energy levels of the Schrodinger equation are 
determined by the zeroes of the Bessel functions Bj with 
integer values j . The TB model gives corresponding results 
at the bottom of the valence band, because the electrons 
are well represented by nearly free electrons. Following the 
ideas presented for the triangles one would expect that the 
energy levels close to the bottom of the conduction band 
could be determined similarly. However, in this case, the 
detailed geometry of the edge (perimeter) has a dominant 
role in determining the energy levels above the Fermi level, 
and the energy spectrum is very sensitive to the number 
of atoms in the circular dot, as shown in figure [2] The 
circular flakes have been obtained by cutting a circle out 
of an inflnite graphene sheet. Note, that the actual edge 
geometry depends not only on the radius but also on the 
site of the center. 

The reason for the size-sensitivity can be traced back 
to the edge states which are present in graphene construc- 
tions with zigzag edges. In circular dots, the perimeter has 
short regions of zigzag segments that are mixed with other 
motifs (especially armchair). This roughness causes edge 
states with different energies, which is in sharp contrast 
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Fig. 2. Geometry of a circular dot cut out from a graphene 
sheet (left) showing electrons density of the state shown with 
an arrow (black is zero density, yellow high density) . The three 
curves (right) show TB-DOS just above the Fermi energy for 
three circular flakes with 3868 (solid line), 3864 (dashed line), 
and 3868 atoms (dotted line), respectively, demonstrating the 
sensitivity of the level structure on edge geometry. 



dots manufactured from semiconductor heterostructures, 
where bound electronic states can exist inside the band 
gap of the semiconductor in question. 
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to the zigzag triangles where all the edge states have ex- 
actly zero energy in the TB model. Figure [2] shows the 
electron density of one such state with density maximas 
at the surface. 

Usually, the shell structure is determined by the overall 
shape of the confining potential in metallic and semicon- 
ductor quantum dots, and the detailed atomic structure 
does not play any role due to the fact that the electron 
wave length is much larger than the interatomic spacing. 
This is not the case in the circular graphene flakes. Al- 
though the wavelenght of the " Dirac electron" is still much 
larger than the interatomic spacing, the tendency for lo- 
calization of electrons close to the zigzag edges destroys 
the simple shell structure, and different circular flakes re- 
sult in qualitatively different electron levels as shown in 
Fig-i 



4 Quantum dots prepared with external 
potential 

So far, we have studied free graphene flakes where the 
electron confinement is determined by the fiake edges. In 
semiconductor heterostructures, quantum dots are usually 
prepared by confining the delocalized conduction electrons 
in a small region with external gates (for a review see [37]). 
It is expected that a similar technique can be applied in 
the future also for the two-dimensional gas of "Dirac elec- 
trons" in graphene. External gates form nearly harmonic 
confinement close to the center of the quantum dot. An- 
other possibility for supported graphene flakes could be 
to modify the atomic structure of the substrate so that 
different regions would comprise different elements, and, 
consequently, cause different interaction with the adsor- 
bate. In this case, the resulting potential well could be 
more of the square-well-type than harmonic. 

The situation is different for an external conflnement 
(infinite graphene sheet) than for a finite fiake with edges. 
Since there is no band gap in graphene, an external poten- 
tial well cannot bind an electron as demonstrated in Fig. 
[3] In addition, this differs considerably from the quantum 




Fig. 3. The upper panel shows an external square well po- 
tential in graphene and schematically the linear energy bands 
in different regions. Electrons cannot be localized by the po- 
tential, because a conduction electron inside the well (dashed 
line) can move out as a valence electron. The two figures in the 
middle display a square-shaped and circular Gaussian poten- 
tial wells in a hexagonal graphene flake. The lower two figures 
illustrate the electron densities of a resonance state above the 
Fermi level (dotted line in the upper panel). 



We have studied the effect of an external confinement 
by using a large hexagonal graphene fiake with 4902 atoms. 
An external potential was added at the center of the fiake. 
We considered three different external potentials: A circu- 
lar well, a square-shaped well, and a smooth Gaussian 
potential with circular symmetry. Surprisingly, the results 
are qualitatively similar irrespective of the type of the at- 
tractive potential. Bound states appear at the bottom of 
the valence band where the electrons act as normal free 
electrons. In the more interesting region close to the Fermi 
level, no bound states can be observed. However, above the 
Fermi level all the potentials result in strong resonances 
with a large enhancement of the wave function amplitude 
within the potential well region. Figure [3] shows the den- 
sities of the wave functions for two such resonances. The 
wave functions do not decay to zero outside the potential 
well but reach a small and uniform amplitude that goes 
all the way to the flake edge (the small amplitude is not 
visible Figure [3]) . 



4 M. Manninen, H. P. Heiskanen, and J. Akola: Electronic shell and supershell structure in graphene flakes 



For a potential barrier, the penetration of a wave func- 
tion inside an apparently forbidden region is often referred 
to as the Klein paradox [34 , which has its origin in the 
Dirac theory of massless fermions. As Figure |3] shows for 
the TB (band structure) model, the wave function pene- 
tration inside the "forbidden region" is a natural conse- 
quence of the missing band gap: An electron that appears 
on the conduction band on one side of the step continues 
as a valence electron on the other side. 



5 Conclusions 

We have studied the possibility of observing electronic 
shell and supershell structure in free graphene flakes. For 
this purpose, we have used a simple tight-binding model 
with one electron per atomic site (pz electrons). Despite 
its simplicity, the TB model describes the key features of 
the graphene band structure close to the Fermi points. 

In large triangular flakes with zigzag edges, the shell 
structure of the "Dirac electrons" in the conduction band 
is the same as for free electrons in a triangular cavity. The 
analytical expression gives the energy levels accurately up 
to ~0.25 eV above the Fermi energy, and the number of 
shells within this region depends on the number of atoms 
in the triangle. A triangle of ca. 40000 atoms (L > 40 nm) 
shows already the first supershell oscillation. 

In the case of circular graphene flakes, the shell struc- 
ture above the Fermi level is dominated by the states that 
are localized close to the zigzag regions of the edges. This 
makes the shell structure very sensitive, not only to the 
radius of the circular flake (number of atoms) but also to 
the location of the center. 

Potential wells which are created on an infinite graphene 
sheet with external potentials (e.g. external gates, inho- 
mogeneous substrate) cannot localize electrons. This is a 
consequence of the missing band gap in the graphene band 
structure. However, such potential wells cause resonance 
states above the Fermi level, which can strongly affect the 
conductance of narrow graphene strips. 
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